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ABSTRACT 


Sixty kHz CW sound was propagated parallel to and near the 
ocean surface adjacent to an array of four thermistors, a sali- 
nometer, a sound velocimeter, a turbulent velocity probe, anda 
wave-height probe in order to investigate the statistics and envi- 
ronmental causes of amplitude modulation. The temporal variations 
of the sound amplitude were studied during four twenty-minute runs 
Over a range of 2 meters at varying depths from near the surface 
to 14 meters during sea-state one conditions in water of 16 meters- 
depth. Analysis indicated mean percentage amplitude modulation of 
approximately 5% with a variance of approximately 1079, Temporal 
correlation times wre about 2 seconds except for a 7 second time 
near the bottom. The power spectral densities of the modulation 
show strong components porecencndnie to predominant ocean wave 
frequencies as well as higher frequencies that may be attributed 
to bubble, temperature and/or salinity inhomogeneities. There is 
evidence that the microstructure dimensions are larger near the 


bottom than in the upper and middle depths. 
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I. INTRODUCTION 


A. ACOUSTIC AMPLITUDE MODULATION IN A RANDOM MEDIUM 

Consider an underwater acoustic wave propagated at a constant 
frequency and a fixed amplitude in isotropically homogeneous, un- 
disturbed water with no scatterers. In the far field of the trans- 
ducer, the wave will arrive at the receiver unchanged except for a 
diminution due to spherical spreading loss and absorbtion. The 
amplitude will not vary about the peak-to-peak value received. 
For low values of absorbtion, this type of acoustic wave is ade- 


quately described by the wave equation: 


or eS: 
WP eS ye 


Near the ocean surface, a wave will encounter random inhomo- 
genieties which cause a varying speed of sound resulting in re- 
fraction. The wave equation must therefore be modified to include 
temporal and spatial changes in the speed of sound. Such changes 
may be expressed as a change in the index of refraction u(x,t), 
where u(x,t) is the ioe of average sound speed to local sound 


speed: Gn 


M(r,t) = 





which may be expressed as: 


aie) a a BORG, G) 





where @ is the root-mean-square variation of the index of re- 
fraction from unity, and n(x ,t) is the normalized variation of U 
from unity. The quantity @ has been shown to be a small number, 
equal to 1.138 x 1074 in laboratory experiments using heated water 
(Stone and Mintzer, 1965) and to 2.55 x 107? in ocean experiments. 
(Whitmarsh, Skudrzyk and Urick, 1957) 

The effect of the random variation of the index of refraction 
on an acoustic wave may be shown by a variation in the amplitude 
of a received signal, noticable as a low frequency modulation of 
the carrier frequency. Such modulation may be expressed as a 
fraction of the maximum undisturbed amplitude of the received 
Signal as shown below: 


VMAX - VMOD 


instantaneous fractional modulation = 
| VMAX 


where 


i 


VMAX the unmodulated voltage value 


VMOD the instantaneous modulated voltage value 


The inhomogenieties which cause the variations in the re- 
fractive index, and ultimately result in amplitude modulation, are 
temporal and spatial variations in temperature and salinity micro- 
structure, random bubble populations, particulate matter, varying 
forms of oceanic life and those inhomogenieties caused by man. 
Fach factor will vary with location and ambient conditions in the 
ocean. The effect on sound transmission may also be a function of 
sound frequency. 

Particle velocity in the ocean has a great effect on the above 


parameters, generally in determining their temporal and spatial 


distribution. There are two main types of oceanic particle velocity, 
that due to convective turbulence and that due to the effect of sur- 
face waves. Convective turbulence is fairly slow and affects the 
distribution of temperature and salinity. Such concepts as tem- 
perature "patches" and lenticular volumes of near unitemperature 

are due to convective turbulence. 

Turbulence arising from the effect of surface wave action is 
generally responsible for random variations in populations of 
scatterers and for small scale variations in temperature and 
Salinity distribution. 

Distribtuions and variations of all parameters affected by 
wave turbulence would show a lessening effect with increasing 
depth due to the 4% decay of wave amplitude. Convective tur- 
bulence could exist at all depths near the surface and its de- 
pendence on depth is unknown. These mechanisms are always present 
near the surface of the ocean and are a principal factor in the 
modulation of sound. They also lead to problems in analysis. For 
example, a stationary temperature patch might show temporal char- 
acteristics due to displacement across a sensor with wave action. 

Statistical analysis of the amplitude modulation occurring in 
inhomogeneous transmission has centered about the coefficient of 
variation, CV, defined by Stone and Mintzer (1961): 


2 2. = 2. 
aes dine = 
7 





where CY = the eoeiricient of Variation 


U 
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the amplitude of the received signal 


oo 
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the ensemble average. 





The same coefficient of variation is defined equivalently by 
Campanella and Favret (1969) as: 


ie 
S(t) — CSlt)7> i 
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where S(t) is the amplitude of the received signal. 

The relationship between the coefficient of variation as de- 
fined by Stone and Mintzer, (CV), and the fractional modulation as 
defined in this thesis 1S shown in Appendix I. 

The relation between (cv)* and the RMS variation of the index 
of refraction @, the range r and the size of the inhomogenieties 
has been variously described, as by Stone and Mintzer (1965) de- 
pending on the (presumed) shape of the spatial correlation function 


of the index of refraction R(p) where: 


K(e) = Cn(FyA(#+eY> 
If R(P) 1S Gaussian, 
2 
kee) = exp (-F 702) 


then the coefficient of variation may be derived for the conditions 


>> 
that a : 
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where K = space wave number of the sound frequency 
Q@ = RMS variation of the index of refraction 
a = radius of the scattering inhomogeniety 
i tiansaueer CO receiver Tange. 











If R(@) is exponential; 
R (eo) = exp (-P/a ) 


then the numerical factor in the expression for (cv)* above would 
be unity. Although an exponential correlation function would not 
satisfy the requirements of non-zero slope at P = O, the experi- 
mental work (Stone and Mintzer, 1961 and 1965, and Campanella and 
Favret, 1969) seems to indicate an exponential correlation 
function. The dilemna can be resolved by uSing a modified ex- 
ponential function to satisfy the Pp = O condition. 

The temporal correlation function for the variation of the 


index of refraction, 


M(ry = (ne) nlt+T)> 


has been compared to the temporal correlation function for the 


variation of sound pressure amplitudes, 
b(t) = CS) StH)? 
and the relation: 


Mird = O(Y) 


which was derived theoretically by Mintzer (1954) and was shown 

to exist experimentally by Campanella and Favret (1969) in a con- 

trolled laboratory environment using convective temperature 

differences to generate a varying index of refraction. 
Experimental work in the ocean environment has included pro- 

pagation at frequencies from 800 Hz (R.M. Kennedy, 1969) and 

60 kHz (Whitmarsh et al., 1957) at distances of 1000 yds 


(Whitmarsh et al.) and twenty-five miles. (Kennedy.) 
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Kennedy found the probability distribution function of the 
logarithm of the amplitude fluctuations to be normally distributed, 
and an exponential temporal correlation function with an average 
decorrelation time to 1/fe in 12 minutes. There was no attempt to 
take meaSurements of the environment during his experiment. 
Whitmarsh, Skudrzyk and Urick found that the fluctuation in amp- 
litude tended to increase with range, and decried a lack of good 


information concerning the microthermal structure of the ocean. 


B. SCOPE AND INTENT OF THE EXPERIMENT 

The purpose of the experiment was to measure and analyze the 
temporal variation in the fractional modulation of an acoustic wave 
Over an Oceanographically defined volume in the ocean. The tem- 
poral characteristics of the environmental factors of temperature, 
salinity, turbulent velocity, sound speed and wave height were 
measured simultaneously with fluctuations in sound amplitude. 

A simple one-transducer and hydrophone arrangement over a 
Short transmission length of two meters was adequate to measure 
Variations in the amplitude of the transmitted signal. The short 
range was desirable because of the need for system rigidity and an 
adequate definition of the environmental factors occurring during 
each twenty-minute run. The length of each run was the minimum 
time necessary to statistically define the low frequency environ- 
mental factors which caused the amplitude modulation. 

The acoustical and Oceanographic analysis at each depth was 
concerned with the statistical and temporal characteristics. The 
results desired were the mean and variance, the probability density 
function, the temporal correlation and the power spectral density 
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The acoustic results will later be cross-correlated with the 
varying environmental factors to determine the relative effect of 


each parameter upon the amplitude modulation. 
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ITI. DESIGN OF THE EXPERIMENT 


A. DETERMINATION OF THE CONFIGURATION OF THE EXPERIMENTAL APPARATUS 
1. Conceptual Design of the Experiment 

In order to measure environmentally caused amplitude modu- 
lation over a clearly defined volume in the ocean, a high frequency 
acoustic carrier wave was directed at a receiving hydrophone 
mounted upon the acoustic axis at a separation of two meters. The 
received signal was processed for the amplitude modulation occurring 
in each twenty-minute run. 

In support of the acoustic investigation, an electro- 
magnetic flowmeter, three thermistors, a Bissett-Berman apparatus 
with a salinometer anda temperature probe and a Ramsey probe with 
a velocimeter and a wave-height probe were mounted about the 
acoustic apparatus to define the environment during each run. 

2. The Acoustic Carrier Wave 

The carrier wave utilized in the experiment was a con- 
tinuous-wave 60 kHz constant amplitude signal. This frequency was 
selected in order to ieeiniee the effect of random bubble popu- 
lations which were anticipated to be greatest near the ocean sur- 
face. 

As may be seen from Figure II-l, the 60 kHz carrier fre- 


quency was near resonance for the most populous bubble sizes near 


me SULtace. The resonance curves were calculated on the basis of 
mae Lormula: {x 
as ane at Re ¥ = Oe. 


13 


Oo 





T-II oeandty 
HLdad SNSHYBA AONANOAYA INVNOSHY Ae dnd 


(zi 19HOTIy) Aouenbary oTysnooy 


OE OOT O06 Og OL O9 OS Ov 
OT 
Zl 
OT 
\ 
V 0) 


\O 





we ve) Yo uy 9 
s Bee we S| [319yM 


(] 9) 
| sntpey etTqqng mos 
sntpey ertqqnd 7109 
sntpey etqqng nol © 
a) 


ST 


N 


es TE fees pean eres eee ance eeay eee ae | 
CO 


© 


Depth (meters) 


14 





where f = resonance frequency (Hz) 


O 
ie radius of the bubble (cm) 

6 
ee = ambient pressure at the bubble depth, equals 10 


, 6 2 
Pes fen at the surface with 10 dynes/cm added 
for each 10 meters of depth 
Pp = density of water (1.0 gram/cm>) 


Y = ratio of specific heats (assumed to be 1.4) 


04 correction for surface tension and thermal 


conductivity equals 1.0296 for 60 kHz 
pe Determination of the Acoustic Ramee 

The most stringent requirements for the distance between 
transducer and hydrophone arose from the desirability that 

(a) the hydrophone be in the far field of the transducer. 

(b) the volume of water under investigation be locally 
defined in terms of the measurable parameters of temperature, 
Salinity, speed of sound and surface wave height. 

(c) the mounting of the hydrophone be stiff enough to 
preclude significant changes of sound amplitude through the sound 
beam. 

(d)imthe path Length be long enough to observe measurable 
amplitude modulation. 

The compromise transducer to hydrophone distance was set 
at two meters. Checking requirement (a), for the available source 
(F-27) using the criterion r = a“[, and assuming an average sound 
speed of 1500 meters a second; 


f2= "60" 107 Hz 


the 250 One meters 


»~ 
lI 


c/f = (1500/6)x 10° 


0.127 meters, where a is the radius of the 


transducer face 


0S) 





then, the far field begins at (approximately) 
r= Aan = OmONGt 27080250 — 0, Olaaemeters 

Thus the far field begins well within the two meter dis- 
tance between the transducer and the hydrophone. 

For requirements (b) and (d) the two meter distance was 
selected to be long enough to yield measurable amplitude modu- 
lation, yet short enough to allow oceanographic measurements on a 
definable volume of water. The data from the runs proved this 
selection to be justified. 

The alignment of the receiving hydrophone with the acoustic 
axis of propagation was an important consideration. A narrowly 
defined beam pattern reduced reflection from the apparatus, but a 
radius of error had to be calculated in order to establish the 
effect of minute movement of the receiver (requirement (c)). 

Pressure amplitude as a function of angle was calculated 
from the formula below, where only the last term is important; 


tea) Zlast-Rr) oie! 
KO sm © 


For one per cent modulation, an angular displacement of 


O ; 
0.45 would have to occur, two meters, the radius of movement 


would be; 
On= 04/05: 
Sin @ = x/2 meters 
<2 (0.0075) 
x = 0.0150 meters 


The radius of movement would have to be 1.50 centimeters, 


more movement than the mounting would allow. As a final test at 
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the scene of the experiment, a diver “"twanged" the hydrophone 
supporting wire with no observable result in amplitude modulation. 
It is estimated that no more than 0.1% modulation could be caused 
by movement of the hydrophone by ocean forces. 

4. Design of the Apparatus 

The apparatus (see Figure IJI-2) as designed was a square 
of height six feet and length two meters from the face of the 
transducer to the receiving hydrophone. 

Three sides of the mount were of two inch diameter steel 
pipe and the fourth was one-eighth inch steel wire stretched by 
means of steel springs to a tension of one hundred and fifty 
pounds. 

The square shape offered: two distinct advantages, the 
shortness of the horizontal members minimized vibratory "tuning 
fork" effects in conjunction with the tension in the wire; it 
also minimized the number of reflection paths from source to 
receiver Acoustic absorbent rubber (SoAB) mounted for eighteen 
inches along the upper and lower members also reduced possible 
reflections to the point that reflections were down 38 db relative 
to the signal when a pulsed 60 kHz signal was propagated. 

An eleven inch steel disc was welded to the vertical pipe 
member directed perpendicularly to the sound axis as a mount for 
the transducer. 

The completed design combined two necessary characteristics, 
rigidity and minimum self-generated hydrodynamic turbulence. The 
end of the apparatus holding the transducer was necessarily 


massive to insure rigidity about the propagating transducer. 
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Hydrodynamic turbulence scattered from the transducer was kept to 
a minimum by having the sound source end of the entire apparatus 
face the mean current. 

The wire served as the mounting for the receiving hydro- 
phone. Since the hydrophone was quite light in weight, the wire 
waS rigid enough, while the small size of the wire precluded any 
significant hydrodynamic turbulence. 

5. The Transducer 

The transducer utilized was the USRD type F-27. The 
transmitting array consists of fifty-five one-inch diameter lead 
metaniobate discs 0.22 inches thick, each of which 1s cemented to 
a one-inch diameter, 4% inch thick tungsten backing plate which 
acts approximately as a circular piston. 

The F-27 transducer was designed for wide range unidi- 
rectional underwater sound projection at a frequency range of one 
kilohertz to forty kilohertz. The USRD recommends a limit of 
forty kilohertz2 since the beam pattern was considered to be too 
sharp at higher frequencies. This characteristic was advantageous 
to the experiment at sixty kilohertz and results were indicative 
of proper transducer sain with the one-half power points 
measured at three degrees off the sound axis, and the side lobes 
at sixty degrees down forty db from the axial sound pressure level. 

6. The Receiving Hy dr ophone 

The experiment required a small, sensitive hydrophone. The 
hydrophone used was the Atlantic Research Corp. type LC-10 with a 
diameter of 0.635 centimeters. For the Sixt yerulonertz carr tex 


wave, the wavelength is 2.50 cm, making the wavelength 3.95 times 


ey 





eivemcGiameter of the hydrophone. The sensitivity of the LC-10 is 
-110 dB re 1 volt per microbar at 60 kHz. The hydrophone was 
mounted vertically along the wire, with the central part of the 


cylinder perpendicular to the sound axis. 


B. PROCESSING ELECTRONICS 
1. Overall System Operation 

The processing electronics were designed to present the 
waveform of the amplitude modulation taking place on the acoustic 
carrier wave over a dynamic voltage range compatible with that of 
the recording equipment. Generally, it was desired to present a 
ten volt peak-to-peak signal at the entry leads of the envelope 
detector, and then process the waveform to view only the amplitude 
modulation. (See Figure II-3.) 

The 60 kHz wave was originated in a GR 1162A frequency 
synthesizer which produced a one volt to four volt peak-to-peak 
waveform accurate and constant to within one millivolt over ex- 
tended periods of operation. 

The carrier waveform was then amplified to the voltage 
necessary to produce a ten volt peak-to-peak waveform at the 
envelope detector by an HP 467A power amplifier set at a fixed 
"times ten" setting. It was found to be advantageous to set the 
power amplifier at a fixed output and vary the output of the 
GR 1162A frequency synthesizer in order to achieve the desired 
voltage value at the envelope detector. 

The power-amplified signal was then transmitted through 
one hundred and fifty feet of waterproof cable to the F-27 trans- 


ducer and transmitted to the LC-10 receiving hydrophone. 
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The Signal as received by the LC-10 hydrophone was amp- 
lified by a fixed +30 db (x 31.6) gain NUS model 2010-030 voltage 
preamplifier in a watertight can mounted on the upper horizontal 
member of the mounting apparatus. The preamplifier required a 
twelve volt DC input, supplied by a Hewlett-Packard 721A DC power 
supply. The amplified signal was transmitted through one hundred 
and fifty feet of watertight cable to the initial electronic pro- 
cessing equipment. 

The initial processing of the modulated waveform consisted 
of bandpass filtering the signal with a center frequency of 60 kHz 
through a pass band of + three hundred Hertz, with -24 db per 
octave attenuation on either side of the pass band. The unit 
iiemhized was the Khron-Hite model 3322 variable filter. In order 
to arrive at a ten volt peak-to-peak waveform at the output of 
the filter, a +20 db internal amplification feature of the filter 
was utilized. The output at this point consisted of a ten volt 
maximum peak-to-peak voltage with a low frequency environmental 
amplitude modulation superimposed upon it, varying to a maximum of 
sixteen per cent of the peak wave form. 

The above oe was processed to detect the modulation 
by the envelope detector (see Figure II-4) described in detail in 
the following section. The salient point is that the ten volt 
fama l was reduced by a factor of 0.6 in the audio transformer of 
the envelope detector. 

The output of the envelope detector was the varying voltage 
induced by the amplitude modulation superimposed upon a small 


amount of DC voltage. 
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CIRCUIT DIAGRAM OF THE ENVELOPE DETECTOR 
Figure II-4 
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DC COMPENSATION VOLTAGE POTENTIOMETER 
Figure II-5 
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The envelope detector output was led into the positive 
Mmoput jack Of a PAR Type 113 differential amplifWer. The dif- 
ferential amplifier was used to eliminate the DC component of the 
amplitude modulation waveform. This was accomplished by simultan- 
eously placing a DC voltage into the negative input jack. The DC 
voltage was supplied by an HP 721A DC power supply set at an 
output of four volts and reduced by means of a 500 ohm potentio- 
meter to that (negative) voltage required to reduce the amplitude 
fMeaulation to a purely alternating Signal. 

The potentiometer was necessary to assure a constant [DC 
input level. The power supply/potentiometer combination was 
tested over a continuous two hour period with less than a one 
millivolt variation during the entire period. The potentiometer 
was configured with an outsized tuning wheel to allow minute (less 
than one millivolt) adjustments of the DC compensating voltage. 
(See Figure II-5.) 

The compensated amplitude modulated signal was amplified 
by a factor of ten in the differential amplifier which made the 
varying voltage too high for the tape recorder, but which was 
subsequently reduced by a potentiometer. The ten-times setting 
was the minimum amplitude setting on the PAR 113. 

The signal was subsequently passed through another Khron- 
Hite 3322 filter set in a four-pole Butterworth low pass configu- 
ration with a high cutoff of ten Hertz and a minimum frequency of 
0.001 Hertz. The attenuation at the high frequency end was -48 db 


per octave. The internal amplification feature was not utilized. 
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The lowpass-filtered amplitude modulation signal was re- 
duced to allow dynamic voltage range compatibility with the tape 
recorder. The reduction was accomplished by the use of a 250 ohm 
potentiometer set to reduce the low-pass filter output by a factor 
Or ten. This reduction had the desirable result of restoring the 
amplitude modulation to its value relative to the six volt output 
of the audio transformer in the envelope detector. 

The signal was recorded on a Sangamo model 3500 fourteen 
track tape recorder simultaneously with the oceanographic mea- 
surements during each run. The dynamic range of the tape recorder 
in the FM mode utilized was 0.01 to 2.50 volts, with all the amp- 
litude modulation observed falling within this range. 

The electronic processing design proved to be quite 
satisfactory except for details in various components which were 
amply compensated for during the course of the experiment. 

aoe tne Envelope Detector 

The envelope detector was adapted from an ocean wave height 
detector and proved to be satisfactory. The detector is essentially 
a double rectification device with an RC filter and a capactive lead 
Heearound potential to discharge DC voltages built up in the sili- 
con diode (1N540) rectifying elements. (See Figure II-4.) 

The time constant for the RC filter at a carrier fre- 
quency of 60 kHz is fifty microseconds and with double rectifi- 
cation there was no appreciable decay between peaks. 

The envelope detector proved sensitive to acoustic modu- 
lation of as little as one twentieth of one per cent in the 


laboratory tanks. 
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Possible improvements to the detector would be the sub- 
stitution of larger capacitors in the filtering element to reduce 


@eessubstantial (1.3 voit) DC component in the output. 
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Lil. “EXPERIMENTAL PROCEDURE 


Ae GENERAL 

The amplitude modulation runs were performed on the twenty- 
first of October 1971 at the NUC Oceanographic Research Tower 
at San Diego. The runs were twenty minutes in length, to allow 
enough time for major low frequency ocean phenomena to be ade- 
quately described and measured. 

The ocean sensors were mounted around the acoustic sensors 
in order to define the volume of ocean over which the amplitude 
modulation took place. (See Figure III-1l.) The sensors were 
selected to measure temporal variations of temperature, salinity, 
turbulent velocity,wave height and the speed of sound during the 
period of acoustic measurement. During each run, the oceanogra- 
phers present also measured and. observed weather and sea surface 
conditions as background for the more precise information pre- 
sented by the ocean sensors. The ocean sensors used are listed 
below with respect to the parameters they measured. 

roa Lenmperature 

a. Three thermistors located about the apparatus obtained 
temporal variation of temperature in the volume of water under 
observation. 
b. One Bisset-Berman temperature sensor measured absolute 

values of temperature near the F-27 transducer over a range of 


14°c to 19%. 
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meeeeoalinity 
One Bisset-Berman salinometer located near the F-27 
transducer measured temporal variation in salinity over a range 
fey to 39.5 parts per thousand. 

Be Particle Velocity 

An electromagnetic flowmeter measured one component, (u) 
or (v), of horizontal velocity and the vertical component (w). The 
flowmeter was located about one-half of the distance between the 
acoustic sensors 1.13 meters downstream from the acoustic axis. 

4. Speed of Sound 

The sound speed was measured as an absolute value with 
temporal variations by a Ramsey Corp. Mark I SVTD probe located 
near the base of the acoustic apparatus. 

5. Wave Height 

The wave height was measured by two sensors, a pressure 
wave height indicator anda Baylor gauge. 

As noted above, not all the ocean sensors were operating 
during all runs. Five runs were made, with runs two through four 
considered suitable for analysis. The general scheme was to per- 
form two runs near the surface, the first run with the electro- 
magnetic flowmeter measuring in the u and v (horizontal) 
directions, and the second with the flowmeter measuring in the u 
and w (or vertical) directions. Three runs followed at deeper 
depths with the flowmeter in the u-w configuration. The first 
run was aborted when the flowmeter malfunctioned, but runs two 
through five proved to be successful. A listing of the runs taken 


1s presented below; all times are local and the date is 21 October 
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1971. Mean water level was 6.7 meters below the concrete deck 
of the NUC tower. The depths are in meters below mean water. 
Roun 1. 
depth: 4.3 meters, 
times: start: 1419 stop: 1439 
comment: The data from this run were not analyzed due to 
a malfunction of the electromagnetic flowmeter. 
eli 2s 
depth: 4.3 meters, 
Memes: start: L530 (stop; 1550 
sensors: turbulence: 1. u (horizontal) velocity 
2. w (vertical) velocity 
3. sound velocity 
4. Both wave heiahts 
temperature: 1. thermistor 7 
2. Bissett-Berman temperature 
Salinity: 1. Bissett-Berman salinometer 
Ue i 
depth: 9.3 meters, 
times: start: 1616stop: 1636 
sensors: turbulence: 1. u velocity 
2. WwW velocity 
Ba Sound Veloce my 
4. both wave heichts 
temperatures: 1. thermistors T, >To; and qT, 
2. Bissett-Berman temperature 


Salinity: Bissett-Berman salinometer 
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Run 4. 
depth: 14.3 meters, 
times: Start: 1045 sop] 170s 
sensors: all sensors operating, see sensor information 
Pie n 3. 
Pan 5. 
depth: 7.3 meters, 
times: start. l7/2omemon mmo 
sensors: all sensors operating, see sensor information 
for Run 3. 
Eee CONFIGURATION OF THE OCEAN SENSORS IN RELATION TO THE ACOUSTIC 

APPARATUS 

Two factors determined the configuration of the ocean sensors 
relative to the acoustic apparatus. First, the volume of water 
around the acoustic axis had to be described adequately for 
temporal variations, and second, the ocean sensors had to be 
placed such that there was little or no interference from down- 
stream turbulence from the sensors themselves. 

Temperature variations and turbulent velocities were considered 
to be the most important parameters affecting the amplitude modu- 
lation. Temperature variations were recorded by the three ther- 
mistors and the Bissett-Berman temperature probe. The Bissett- 
Berman probe was placed centrally between the acoustic apparatus 
and the mount for the flowmeter due to its size. Two of the ther- 
mistors were placed above and below the receiving hydrophone in the 
vertical plane of the acoustic axis, while the third was placed 
along the horizontal plane of the axis midway between the trans- 


ducer and the hydrophone. (See Figure III-2.) 
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The electromagnetic flowmeter was placed so that it measured 
the turbulence occurring 113 centimeters from the center of the 
ee@mstic axis. 

All other sensors were integral parts of large instruments 
and were placed between the acoustic apparatus and the mount for 
the flowmeter to minimize the effect of their wake turbulence. 

aur e III-3 shows the distances between various sensors. The 
mean current was as shown in Figure III-l, with the acoustic 


apparatus up-stream of all other equipment. 


C. DETERMINATION OF RUN DEPTHS 

The mean water depth at the NUC tower was 16.4 meters (54 feet) 
during the time of the runs. Four depths were selected in order to 
observe near-surface to near-bottom amplitude modulation effects. 
The first two runs were at an estimated water depth of 4.3 meters, 
the third at an estimated 9.3 meters, the fourth (and deepest) at 
14.3 meters, and the last run WAS Ab je Meterce 

These depths were chosen in order to cover the entire water 
column and to minimize interference from the underwater beam 
Structure of the tower. The shallowest depth was far enough from 
the surface to avoid undue stresses on the apparatus, while the 
deepest depth was far enough away from the bottom to avoid inad- 


Vem@tamtc contact. 


D. THE NUC TOWER 
The NUC tower is located approximately one mile off Mission 
Beach, California in 16 meters of water. The apparatus as de- 


picted in Figure III-2 was mounted on the west side of the tower. 
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The electronics were installed on the main deck, while the appa- 
ratus was handled from the second deck immediately below. (See 
Figure III-4.) 
Environmental factors were observed at the tower while the 

runs were in progress. The factors noted were: 

1. swell: from the northwest, two feet from crest to 
trough, with an estimated period of four to five seconds. 

2. wind: 310°, 300°, 320° from 9 knots to 7.5 knots during 
the course of the experiment 

3. wind waves six inches from crest to trough 


lk. some bubbles visible on the surface 


Ee UNUSUAL EFFECTS 

Unusual effects were those effects arising from considerations 
Other than physical parameters and naturally occurring bubble 
populations. All the unusual effects noted were due to animal and 
plant life in and around the apparatus during each run. 

The area around the tower supported a large population of 
small fish ranging up to eight inches in length. It was noted 
that fish were attracted to the apparatus when it was lowered 
initially, but resumed their normal habits after a short period 
of time. Fish could not affect the wire mounted hydrophone, nor 
were any effects noted, since fish eee through the sound beam 
would affect the modulation for a minimal amount of time. 

Piant life could attect the equipment by criftting imtoo it. 
There was a small amount of kelp in the area at all depths, but 


no acoustic effects noted during the runs. 
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Boat traffic around the tower had some effect, and could have 
been sources of turbulent velocity and bubbles which could noti- 
mil affect the amplitude modulation of the acoustic signal. This 
effect was quite apparent during the Run 5, when there was a boat 


idling in the area. 


F. ORGANIZATION OF INDIVIDUAL RUNS 

All information from each Sensor was recorded on Sangamo 
Model 3500 tape recorder. Between individual runs, all sensors 
were monitored for proper output. After a depth change, a few 
minutes were allowed for all instruments to be checked for proper 
Operation. When all instruments were checked and operating, the 
tape recorder was started. The simultaneous start and stop times 
controlled by the tape recorder insured an accurate run period 
menmail Sensors. 

The recorder channel allocation was as follows: 


Channel 1.: Ramsey probe temperature 


Channel 2.: Ramsey probe velocimeter 
Channel -: Bissett-Berman salinometer 


Channel 


3 

Channel 4.: Turbulent velocity (v) 
Ds? Purbulent veloot cya 
6 


Channel 6.: Baylor wave height probe 
Channel 7.: Amplitude modulation 
Channel 8.: Pressure wave height probe 
Channel 9.: Thermistor T, 


Channel 10.: Bissett-Berman temperature 


Channel 11.: Thermistor 7, 
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Channel 12.: Thermistor T. 


Channel 13.: Phase modulation 


Channel 14.: Voice 
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TV. “ANAEYSTSsOm ACOUSTIC RESULTS 


fee LNITIAL DATA REDUCTION 
1. General 

The recorded analog data consisted of four twenty minute 
runs at varying depths. The dynamic range of the recorded voltage 
was 0.01 to 1.10 volts. The data reduction scheme involved three 
steps; transfer of the analog data from magnetic tape to a Brush 
Mark 200 strip-chart recorder, conversion of the printed analog 
fata tO digitized data and computer analysis of the digitized data. 

mee Magnetic Tape to Strip-=Chart Record 

The tape recorded data was originally collected at a tape 
speed of 1 7/8 inches per second. When replayed for the Brush 
Mark 200 eight channel strip-chart recorder, the data was played 
at 7% inches per second, making the strip-chart data for each run 
one fourth as longo as the real-time runs. This had the effect of 
compressing a twenty minute run to five minutes of strip-chart 
readout and increasing the apparent frequencies of modulation by a 
jeer Of four. At MP ae recorder speed of five milli- 
meters per second, each printed run was presented on 1.5 meters 
Stee nart. 

The strip-chart recorder was set to present the data at 
50 millivolts per millimeter, or 1.27 volts per inch. 

The other seven strip-chart channels were used to visually 
compare some of the ocean Sensor outputs to the output from the 


amplitude modulation data. The other seven strip-chart channels 


By 





were allocated as follows: 
Channel 1.: Salinometer 
Cha nels 2 2 ev veloc ie, 
Channel 


-: w velocity 


Channel 


S) 

Channel 4.: Amplitude modulation 
5.: Baylor wave height 
6 


= Thermistor Lf 


3 


Channel 7.: Bissett-Berman temperature 


Channel 


Channel 8.:; Velocimeter 

Jt was then possible to make preliminary visual comparison 
of the waveforms of the recorded outputs in order to observe any 
Similarities in the waveforms. (See Figures IV-1 through Iv-4.) 

prot ip -Chart kecord te Digital Tape 

The strip-chart records of each of the four runs to be 
analyzed were taken to the Fleet Numerical Weather Central facility 
ee oint Pinos, California to be recorded as digitized data on 
magnetic tape. The data were digitized directly from the strip 
Poet ©TeCCOLGiIng om the Pacruity,, Ss tracing Gigitazer. The digi- 
tizer recorded data points in increments of hundredths of an inch 
in the y-direction for each one one-hundredth of an inch of travel 
in the x-direction, and transferred the x and y values to a nine- 
Ueeck magnetic tape. 

Thesea si tizing Wacmamitiatedeateunesstami mol each run, 
visible on the strip-chart as a vertical line produced when the 
tape recorder run was commenced. The zero amplitude value in the 
y-direction and the zero time value in the x-direction were taken 


at this line. There were data excursions in the positive and 
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EIGHT CHANNEL BRUSH RECORDER READOUT FOR RUN 2, DEPTH: 
4.3 meters. top to bottom: salinometer, v-turbulence, 
w-turbulence, amplitude modulation, wave height, ther- 
mistor 3, Bissett-Berman temperature and velocimeter 
readouts. 


Figure IV-l 





EIGHT CHANNEL BRUSH RECORDER READOUT FOR RUN 5 DEPTH: 
7.3 meters. For channel identification, see Fig. IV-l. 


Figure I[V-2 








EIGHT CHANNEL BRUSH RECORDER READOUT FOR RUN 3 DEPTH: 
9.3 meters. For channel identification, see Fig. IV-l. 


Figure IV-3 





EIGHT CHANNEL BRUSH RECORDER READOUT FOR RUN 4 DEPTH: 
14.3 meters. For channel identification see Fig. IV-l. 


Figure Iv-4 





negative directions, measured in one-hundredths. The conversion 
factors were 0.0127 volts per data point in the y-direction, and 
0.01265 seconds (real time) per data point in the x-direction. 
4. Digitized Data_from Magnetic Tape to IBM Data Cards 

The conversion of the digitized data from magnetic tape 
to data cards for analysis on the IBM 360 digital computer was 
accomplished in the Fleet Numerical Weather Central's CDC 5000 
computer. 

The procedure utilized Program I, and converted at least 
5600 data points for each run, corresponding to nineteen minutes 
of real data collection time. This decrease in data time was due 
to an uncertainty of stopping points on the strip-chart readouts, 


and allowed an identical number of data points for each run. 


beeen NALYSIS OF ACOUSTIC DATA 

The data from each run were analyzed on the Naval Postgraduate 
School's IBM 360 digital Series 

The analyses were performed on time-varying amounts of 
fractional amplitude modulation. The raw data as presented on 
data cards was in varying values expressed as positive and 
negative hundredths of an inch. The conversion was common to all 
programs and is described below: 

a. Each run was scanned for the maximum positive value. 

b. The maximum positive value was subtracted from all other 
values and converted to 0.0. This level corresponded to the 
unmodulated amplitude. All data values were now negative values, 


all less than zero except for the maximum value. 
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em bach Value was multiplied by a conversion factor consisting 
of (-1.27) which converted the data points to positive voltage 
values, and (1/6) which made each point a fractional value of the 
maximum amplitude of six volts. The resultant data points could 


now be expressed as fractional modulation: 


f= (unmodulated voltage) - (modulated voltage) 


(unmodulated voltage) 


The converted raw data was then analyzed for the mean, 
variance, temporal correlation and power spectral density of each 
run in Program II. Program JII computed a histogram of the amount 
of fractional modulation versus the number of occurrences of each 
value. An independent calculation of the mean, mean Square and 
variance was made in Program IV. (See Figures IV-5 through IV-7.) 
The results of the various analyses are presented in the following 


Section. 
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DATA ANALYSIS SCHEME: PROGRAM IT 


RAW DATA ON IBM CARDS 





FIND MAXIMUM VALUE (XMAX) VALUE OF UNMODULATED 
SET EQUAL TO ZERO AM nD ee Oe 


emeBTKACT ALL VALUES FROM | (NEGATIVE) VALUES OF XMAX 
MAX VALUE (X, -XMAX) PESoeeA Che HODULATE PD Value 


= 


CONVERTS ALL VALUES TO 
ElerrPLy VALUES BY (-1.27) POSITIVE VOLTAGE DIFFERENCE 






a freon 


DIVIDE ALL VALUES BY (6.0) 


a ae 


CONVERTS ALL VALUES TO PER- 


hee MODULATION 


PAstnvs LOUNT  ALIioC eV AZURE S Maybe EXPRESSED Ao: 
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Cn ee 
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FIND VAR TANCE 
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Figure IV-5 





DATA ANALYSIS SCHEME: PROGRAM ITI 





RAW DATA ON IBM CARDS 


CONVERSIONSOF ALL DATA VALUES 


TO FRACTIONAL MODULATION VALUES 






ESTABLISH "BINS'' FROM 0.0 TO 0.170 

IN STEPS OF 0.001, YIELDING 171 ''BINS' 
(0.160 IS THE MAXIMUM VALUE OF 
FRACTIONAL MODULATION) 





SCAN EACH VALUE OF FRACTIONAL MODULATION 
AND PLACE IN THE APPROPRIATE BIN 





COUNT THE NUMBER OF OCCURENCES IN EACH BEN 


Figure IvV-6 
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RAW DATA ON IBM CARDS DATA ANALYSIS SCHEME: PROGRAM IV 


SeOnveRSION TO F 


as > PRINT MEAN VALUE 


CALCULATE ( F, ) 


PRINT MEAN SQUARE VALUE 





=r 
CALCULATE | A 


sz 


Is |. CE) V} 


Figure IV-7 





PRINT SIGMA SQUARE 
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V. ACOUSTIC RESULTS 


fee LHE STRIP=CHARIT READOUTS 

Figures V-1l1 through V-4 are representative 168-second portions 
of each twenty-minute run. The variation in voltage from peak-to- 
peak 1s approximately One volt, and the apparent frequencies of 
variation are four times the actual frequencies occurring during 
each run. 

From visual inspection of the strip-charts, it may be seen 
that the modulations occurring during the two shallow runs (4.3 
and 7.3 meters) were more active in frequency than during the two 
@eemer runs (9.3 and 14.3 m.). Lt must be noted that ere og meter 
run was conducted after the other runs, and that during most of the 
period a fifty-foot motor launch was in the area of the tower with 
with propellors turning. The launch served as a bubble injection 
device as well as a turbulence generator. The run was not invali- 
dated because the modulation effect of bubbles artificially gen- 
erated during this run could be compared with that of the ambient 


medium of the previous runs. 


B. MEANS AND VARIANCES 
The mean and variance of the amplitude modulation for each run 
was determined by Computer Run IV. The formulation was based on 


the fractional modulation F, and was calculated as shown: 


= = | YMAX_- VMOD 
VMAX 
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BRUSH RECORDER READOUT OF AMPLITUDE MODULATION RUN TWO DEPTH: 
4.3 meters. vertical lines indicate four seconds in real time, 
horizontal lines indicate 0.25 volts. 


Figure V-l 
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BRUSH RECORDER READOUT OF AMPLITUDE MODULATION RUN 5 DEPTH: 
7.3 meters. see caption for Fig. V-l. 


Figure V-2 
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BRUSH RECORDER READOUT OF AMPLITUDE MODULATION RUN 3 DEPTH: 
9.3 meters. See caption for Fig. V-l. 


Figure V-3 





BRUSH RECORDER READOUT FOR AMPLITUDE MODULATION RUN 4 DEPTH: 
14.3 meters. see caption for Fig. V-1l. 


Figure V-4 
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which resulted in the time series F(t); then, 


Cece = et 


al 


where n = 5600 


a CEA > = Zr 


Yn 


and oe KK FUE)? - gece $ 


The results obtained from each run were compared with the 
results of Computer Run II, which determined the mean, temporal 
correlation and the power spectral density for each run. Effecti- 
vely, for Program II the mean, (F(t)), was determined as in 
Program IV so that comparison was superflous. 

ihe Variance, a”, can be formulated in three ways; first as 


shown above, secondly as the zero-time temporal] correlation value, 


TA. 
ae : 2. 
— an Cue ( 
6 Sd Md ie ae OF CEV Fler rate 
{< a Vo 


and thirdly as the value of the area under the curve for the power 


spectral density: 
z 
6 fE(Fid 


It should be noticed that at all depths, the three values of 
the variance are quite close. The slightly different values of 
a“ from the calculation for C(O) and from the simple sorting cal- 
eulation (Program IV) are solely due to round off error for the 


matter -calculation. 


Sill 





The value of o* derived from the integral of the power spectral 
density has the largest numerical value. This is the result of 
limiting the higher frequencies from the calculation of the Fourier 
fers Dand, resulting in a rounded "window". This rounding off 
resulted in a very Slightly larger area under the power spectral 


density curve than would otherwise occur. 


The results for each run are shown below: 


Run 2: 4.3 meters 
Program IV: (F) = 0.0526 
2 
Chae o.004 1 
a” = 0.0013 
Program eG) = 0) O52 16a 
2 
i= Cl Oa OF OOmro 1 
pas 
0 = SG(r) GF = 0.0014469 
loum 5: 7.3 meters 
Program IV: (F) = 0.0460 
2 
Cem) a 40.063) 
og? = OO LO 
Program II: (F) = 0.04604 
2 | 
O”~ = C(O) = 0.00096 
ic 
o° =JS E(F)dF = 0.0009756 
Run 3: 9.3 meters 
Program IV: (FE) = MOE Os 
(F“) = 0.0016 
of” = ~2OO7 
Program II: (F) = 0.03026 
Zan = is 
6" ="6(0), = O-000G, 
2 Q 
o” = SEe(F)dF = 0.0007378 





Run 4: 14.3 meters 


Program IV: (F) = 0.0470 
(r~) = 0.0040 
o” = 0.0018 
Program II: (F) = 0.0472 
o” = C(O) = 0.00173 
o” = SE(F)dF = 0.0019268 


The probability density of fractional modulation values for 
each run is shown on the histograms. (Figures V-5 through V-7.) 
Although the mean values calculated above do not differ in any 
orderly manner with depths, the histograms (as well as Figures 
V-1,2,3 and 4) show that for the shallow depths there are greater 
numbers of short duration, high modulation values, whereas for the 
greater depths, the distribution of fractional modulation centers 
heavily about the lower values. The bimodal P.D.F. at the greatest 
depth (14.3 meters) shows that onset of infrequent occurrences of 
very high modulation values. This is the statistical representation 
of the type of modulation observed on Figure V-4, where the modu- 
lation assumes a high value and varies about that value for 
relatively long periods. The data thereby suggest that as depth 
1s increased the magnitude of the medium inhomogenieties increase 
and the size (temporal extent) of the inhomogeneous volumes becomes 


greater as well. 


See leMPORAL CORRELATION 
The temporal correlations (time autocorrelations) for each run 


are shown on Figures V-8 through V-ll. The times to 50% 
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FRACTIONAL MODULATION 
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14.3 meters 


Fractional Modulation 
Figure V-7 
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NORMALTZED TEMPORAL CORRELATION OF RUN 2 
DEPTH: 4.3 meters 


Figure v-8 
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NORMALIZED TEMPORAL CORRELATION OF RUN 5 
DEPTH: 7.3 meters 


Figure V-9 
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Fractional Correlation 
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NORMALIZED TEMPORAL CORRELATION OF RUN 3 
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Figure V-10 
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NORMALIZED TEMPORAL CORRELATION OF RUN 4 
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Figure V-11 
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Seemrelation and l/e correlation are shown below: 


Time to 50%: Depth (meters) Time (seconds 


L3 Piss s' 
Te eS 
Se 5, 223 
14.3 6.9 
Time to l/e: ues Bae 
03 1.65 
ons 2S 
14.3 Te) 


Tt can be seen from the figures that there 1S a rapid initial 
decorrelation down to low as 10% for depths of 4.3 and 9.3 meters, 
and a slower decorrelation for the 14.3 meter depth. All of these 
runs show a fluctuation between steadily decreasing minima and 
maxima with a delayed zero crossing. Whitmarsh, Skudrzyk and 
Urick (1957) hypothesized that amen fluctuation was an indication 
of the size of the inhomogeneous "patches" around the receiver. 

If Whitmarsh, Skudrzyk and Urick's hypothesis 1s correct, then the 
apparent "patch" sizes show an increase with depth. If this is the 
case, then small-Scale turbulence may tend to dominate at all depths 
near the surface. 

The temporal correlation of the 7.3 meter depth is of interest 
Since it 1s known that an artificially derived bubble and turbulence 
frend was introduced into the volume under consideration. The 7.34 
meter depth shows the most rapid decorrelation to zero of all ee 


and twice as frequent fluctuations about zero as for the undisturbed 


Gk 





a a 





water at other depths. This suggests that the ambient "patches" 

of water were chopped into smaller patches by the boat's propellor. 
Another hypothesis might be that the inhomogeneous volumes 

are more or less "frozen" and that the spatial displacement of 

these volumes about the acoustic axis may cause the appearance 

of temporal variation. Corroboration of these hypothetical 

explanations awaits the complete analysis of the data gathered by 


the oceanographic sensors. 


D. POWER SPECTRAL DENSITY 

The power spectral densities (PSD) of each run are plotted 
against frequency on Figures V-12 through V-15. The power spectral 
densities show the highest values at those frequencies (0.08 to 
O.12 Hz) which appear to correspond with the observed ocean swell 
frequencies prevailing during each run. 

The PSD's show a diminishino value both with increasing fre- 
quency and with increasing depth. The two shallowest depths 
Ga3 and 7.3 meters) show a similarity in the values of PSD from 
wee .O Hz, except for a high "DC" value in the plot for 4:3 
meters, resulting in a higher contribution to the integrated value 
oa? SD (0~) at that depth. A similar likeness is apparent between 
the two greater depths (9.3 and 14.3 hecouee For frequencies 
greater than 0.1 Hz, the PSD at the two shallower depths appear 
to be consistently four to five db greater than for the deeper 
depths. 

The frequencies of relative maxima along each of the PSD curves 
show similarities at each end of the spectrum. At the low fre- 


quency end, the ocean swell frequencies are more pronounced at 
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the shallower depths of ie 2 anG@e7 33 Meterc. ae inece no«ima are 
Pioeher tor the’ 7.3 meter run than for the 4.3 meter run, which 
may again show the effect of the boat in the area injecting 
bepbiecs and turbulence into the medium. At 9.3 and 14.3 meters, 
these low frequency inaxima tend to dissappear into the overall 
curve, showing a decrease of the component of the amplitude modu- 
Mee1on due primarily to the turbulent effect of surface waves. 
The amplitude modulation at the higher frequencies (0.2 to 
1.0 Hz) could be due to spatial distributions of temperature, 
Salinity or bubbles which are seen as temporal frequencies due to 


momecment across the acoustic axis. 


Eee CONCLUSIONS 

The conclusions presented here await the additional infor- 
mation that the data from the ocean sensors will provide. This 
information will be available in the theses of LCDR Duchock and 
_T's Seymour and Bordy, to be ce at the Naval Postgraduate 
School in March, 1972. The present conclusions are based upon 
the acoustic data alone and are listed below: 

1. Amplitude modulation near the ocean surface ranged in 
value from zero to sixteen per cent. The maximum amplitudes were 
not noticably restricted to any one depth, but as depth increased, 
the distribution of modulation values tended to center about the 
lower values. 

2, The temperal correlation eumves for the #.3, 924, ancdelos 
meter runs indicate that the size of the volumes of inhomogenicties 


("cellis'') increases as the depth 1S inerecased. The correlation 
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curve for the 7.3 meter run shows a significant decrease in the 
Size of the volumes, a direct result of having a motor launch in 
the area injecting turbulence and bubbles into the volume of water 
under investigation. Thus it may be seen that in the ambient 
ocean, cells of microstructure are larger at greater depth. Near 
surface cells would tend to be reduced in size by the action of 
turbulence arising from the effect of surface waves. 

3. A comparison of the variances of the two shallow runs, 
fae. 5 and 7.3 meters, shows that there is a small) difference 
between the two. (0 (4.3) = 0.0013 and G-(7 a = 0.0010.) Com- 
paring the temporal correlations for the two runs it maybe seen 
mee the decorrelation time for the 7.3 meter run is One-half that 
Semthe 4.3 meter run, with more numerous fluctuations about the 
@emoeaxis. Such data suggests a reduction in the size of the cells 
of inhomogeniety. 

Mintzer's formula for the variance of acoustic amplitude 


modulation, derived for an exponential correlation function (see 


Appendix TI) 2 1 9 2 
(cy) qr Ky ar 


shows the dependence of variance on the wave number, k, range, r, 
@meeRMS variation of the index of refraction, @, and the size of 


the inhomogenieties, a. In this experiment, wave number and range 


° 
a 


were constant; since neither GO nor (cv)* changed appreciably, it 
can be deduced that there was a simultaneous increase in @ (with 
Pae reduction of a) for the 7.3 meter run. 

4, The parameters which could cause a change in the RMS 
Yarmation of the index of refraction, &, Were wenperaiimc and 


Scmmiaty microstructure, and bubbles. 
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The effect of resonant bubbles on the change in & would 
be very great. Since bubble populations decrease with depth, this 
might account for the implied decrease in @ with depth. When in 
the presence of a motor launch, the size of the inhomogeneous 
volumes were made smaller, there was a concurrent increase in the 
bubble population and in the value of @. The changes and vari- 
ations in the other parameters, chiefly temperature and salinity, 
should be accounted for in subSequent analyses, and a quantitative 
estimate of the relative effect of bubbles on amplitude modulation 
may be possible. 

5. The power spectral densities indicate that most of the 
modulation takes place at those frequencies corresponding to the 
frequencies of the surface waves. This effect was noticable down 
to the greatest depth (14.3 meters) although the prominemce of 
the PSD's at the surface-wave frequencies tended to disappear with 
greater depth. 

6. In summary, the major causes of acoustic amplitude modu- 
lation at 60 kliz in this near-surface volume of ocean could well 
have been the effects of mobile, random bubble populations acted 
upon by particle velocities generated by the action of surface 


Waves. 


F. SUGGESTIONS FOR FURTHER EXPERIMENTS 

The basic questions left incompletely answered by this 
investigation are the relative effect of bubbles versus the other 
causes Of amplitude modulation near the surface, and the spatial 
effect of the inhomogenieties upon the transmission of sound. 
Later analysis of the oceanographic data will help to answer this 


question. 
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It 1S Suggested that further experiments of this kind would be 
valuable. In order to investigate the spatial effects of the 
inhomogenieties, the basic apparatus of this experiment could be 
modified as shown in Figure V-16 in order to attempt some cross- 
correlation of the amplitude modulation received by each of the 
hydrophones. In order to more completely investigate the question 
of transport versus the stationary oscillation of inhomogenieties, 
more thermistors could be mounted as shown on Figure V-16. 

The question cf bubbles and bubble populations in an area 
might be answered either by analysis of all other parameters and 
theix effect, or by deliberately generating bubbles in the volume 


mmEvAter to note their effect on ambient conditions. 
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APPENDIX 
THE RELATION BETWEEN F AND CV 


The relation between F, the fractional modulation which is 
used in the development of the experiment, and CV, which Stone 
and Mintzer call the "coefficient of variation" will be developed. 

Fractional modulation is defined in terms of the voltages as 


shown below: 
— Ee Y | 
VMAX 


where V = modulated voltage 


I 


VMAX unmodulated voltage 


The coefficient of modulation is defined by Stone and Mintzer 


in terms of pressure amplitude, P, as shown below: 
a 2 wy 
(Cv) = (<pty = SP2 3 
ie 
<P> 


Then i1f one assumes that P = KV, or that a constant factor 
relates the two values, as it does in the equipment used for these 


experiments, then: 


> At a 
h iv e a CeVP- } ae Cyt? on eve 


oe 


LN . 





Fractional modulation was used in this experiment because 
small fluctuations were anticipated. <A development between the 

; ; aa 
variance of the fractional modulation and the term (CV) is shown 


below: Be vee y | = (1 Da \ 
VAX 





Visa | 


d 





IKre>-<F? | 
ae ‘ile Yul 7 eae : 
ia awe — (<1 - a \ 


3 


ra 


’ "2 
- 24e » 4 c( me) > . <1> + 4 ro ~ (iuew? 


<< * ap ) ia any 


or, Since VMAX rena weOuecma nts 
4 
2%, oo ye Kv? 
4 aw 
(Tern) 
. lve 
moet aolied by the factor Bes the variance may now be 
Sve 
expressed: 2 
CNM iis : 
j cv > : etree nant amen ann ann rare 
Nee 
(\inax)™ <v ; 


Where the second term 1s a numerical value multiplied by 


Mintzer and Stones' expression for the coefficient of variation. 


The value multiplied against the coefficient of variation is 


<v? 


Te 


(YJ AX s 


which is numerically obtainable as follows: 
CED>s oe 
a <p a ee, 


mor all runs, the value of ¥ was 6 Volts, 





MAX 
iC) ee Ce 
(Vesey) * 3b 


WS 





iecett 
The relation between ae and (CV) will be calculated for the 
4.3 meter run: 


2 ce 
ve (ie Oosaed Codd og 
(yaar)? 36 : or oe 





» 


and in the 7.3 meter run was: 


2. or Sin Meee 
év> (1~ 0.0460!) . (0.9840. 9,4268 
a ea ae ee ee Loe Fy = 
(vice) cs SG 


It 1s noted that the conversion factor does not significantly 
change with depth, and therefore conclusions about the relative 
zZ : : 
value of (CV) at two different depths are applicable as well as 


; 2 
the relative value of 0 at these depths. 
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